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FUNCTIONAL CENTRAL LIMIT THEOREM FOR SUBGRAPH COUNTING 

PROCESSES 

TAKASHI OWADA 


Abstract. The objective of this study is to investigate the limiting behavior of a subgraph counting 
process. The subgraph counting process we consider counts the number of subgraphs having a 
specific shape that exist outside an expanding ball as the sample size increases. As underlying 
laws, we consider distributions with either a regularly varying tail or an exponentially decaying 
tail. In both cases, the nature of the resulting functional central limit theorem differs according 
to the speed at which the ball expands. More specihcally, the normalizations in the central limit 
theorems and the properties of the limiting Gaussian processes are all determined by whether or 
not an expanding ball covers a region - called a weak core - in which the random points are highly 
densely scattered and form a giant geometric graph. 


1. Introduction 

The history of random geometric graphs started with Gilbert’s 1961 study m) and, since then, 
it has received much attention both in theory and applications. More formally, given a finite set 
A C M'’* and a real number r > 0, the geometric graph G{X, r) is defined as an undirected graph 
with vertex set A and edges [x,y] for all pairs x,y £ X for which ||x — y|| < r. The theory of 
geometric graphs has been applied mainly in large communication network analysis, in which the 
connectivity of network agents strongly depends on the distance between them; see m, [26], and 
Chapter 3 of HZ). On the purely theoretical side of random geometric graphs, the monograph [2T] 
is probably the best known resource. It covers a wide range of topics, such as the asymptotics of 
the number of subgraphs with a specihc shape, the vertex degree, the clique number, the formation 
of a giant component, etc. From among these interesting subjects, the present study focuses on 
constructing the functional central limit theorem (FCLT) for the number of subgraphs isomorphic 
to a predefined connected graph T of hnite vertices. 

A typical setup in |21] is as follows. Let A„ be a set of random points on Typically, this will 
be either an i.i.d. random sample of n points from /, or an inhomogeneous Poisson point process 
with intensity nf, where / is a probability density. We assume that the threshold radius depends 
on n and decreases to 0 as n ^ oo, but we do not impose any restrictive assumptions on / except 
for boundedness. Then, the asymptotic behavior of the subgraph counts given by 

(1-1) := X = r} , 

ycXn 

(= denotes graph isomorphism, and T is a fixed connected graph) splits into three different regimes. 
First, if nrjj 0, called the subcritical or sparse regime, the distribution of subgraphs isomorphic to 
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r is sparse, and these subgraphs are mostly observed as isolated components. If nr'^ —€ (0, oo), 
called the critical or thermodynamic regime, for which decreases to 0 at a slower rate than the 
subcritical regime, many of the isolated subgraphs in G{Xn,rn) become connected to one another. 
Finally, if nr^ —>■ oo (the supercritical regime), the subgraphs are very highly connected and create 
a large component. 

Historically, the research on the limiting behavior of subgraph counts of the type (jl.ljl dates 
back to the studies of m, m, and [27], in all of which mainly the subcritical regime was treated. 
Furthermore, [7| adopted an approach based on the martingale CLT for [/-statistics and proved a 
CLT under various conditions on / and r„. Relying on the so-called Stein-Chen method, a set of 
extensive results for all three regimes was nicely summarized in Chapter 3 of m- Recently, as a 
higher-dimensional analogue of a random geometric graph, there has been growing interest in the 
asymptotics of the so-called random Cech complex. See, for example, [18], [T9], and [28], while [TO] 
provides an elegant review of that direction. 

Somewhat parallel to but more important for the study on the geometric features of 

extreme sample clouds, is an alternative that we explore in this paper. To set this up, we introduce 
a growing sequence Rn —>■ oo and a threshold radius t > 0. The following quantity, Gn{t) counts 
the number of subgraphs in G{Xn,t) isomorphic to F that exist outside a centered ball in with 
radius Rn- 

(1.2) Gn{t) := l{G(T,t) ^ F} X l{m{y) > Rn} , 

where m(xi,... ,Xk) = mini<j<fc ||xj||, Xi € and || • || is the usual Euclidean norm. 

From the viewpoint of extreme value theory (EVT), it is important to investigate limit theorems 
for Gn{t)- Indeed, over the last decade or so there have been numerous papers treating geometric 
descriptions of multivariate extremes, among them a, 0, and [6]. In particular, Poisson limits of 
point processes possessing a U-statistic structure were investigated by [12] and [2l] , the latter also 
treating a number of examples in stochastic geometry. The main references for EVT are [14] . [22], 
and [13] . 

The asymptotic behavior of (jl.2|) has been partially explored in [20], where a growing sequence 
Rn is taken in such a way that has Poisson limits as n —)■ cx). The main contribution in [20] 
is the discovery of a certain layered structure consisting of a collection of “rings” around the origin 
with each ring containing extreme random points which exhibit different geometric and topological 
behavior. The object of the current study is to develop a fuller description of this ring-like structure, 
at least in a geometric graph model, by establishing a variety of FCLTs which describe geometric 
graph formation between the rings. 

By construction, the subgraph counts (|1.2h can be viewed as generating a stochastic process in 
the parameter t > 0, while a process-level extension in (II.Ih is much less obvious. Then, while 
m captures the dynamic evolution of geometric graphs as t varies, dni) only describes the static 
geometry. Thus, the limits in the FCLT for (II.2h are intrinsically Gaussian processes, rather than 
one-dimensional Gaussian distributions. 

One of the main results of this paper is that the limiting Gaussian processes can be classihed 
into three distinct categories, according to how rapidly Rn grows. The most important condition 
for this classification is whether or not a ball centered at the origin with radius Rn, denoted by 
B{0,Rn), asymptotically covers a weak core. Weak cores are balls, centered at the origin with 
growing radii as n increases, in which the random points are densely scattered and form a highly 
connected geometric graph. This notion, along with the related notion of a core, play a crucial 
role for the classihcation of the limiting Gaussian processes. Indeed, if B{0,Rn) grows so that it 
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asymptotically covers a weak core, then the geometric graph outside B{0, Rn) is “sparse” with many 
small disconnected components. In this case, the limit is denoted as the difference between two 
time-changed Brownian motions. In contrast, if B{0,Rn) is asymptotically covered by a weak core, 
the geometric graph in the area between the outside of B{0, Rn) and inside of a weak core becomes 
“dense”, and, accordingly, the limit becomes a degenerate Gaussian process with deterministic 
sample paths. Finally if B{0,Rn) coincides with a weak core, then the limiting Gaussian process 
possesses more complicated structure and are even non-self-similar. 

We want to emphasize that the nature of the FCLT depends not only on the growth rate of Rn 
but also the tail property of /. This is in complete contrast to because, as seen in Chapter 

3 of [21], the proper normalization, limiting Gaussian distribution, etc. of the CLT are all robust 
to whether / has a heavy or a light tail. In this paper, we particularly deal with the distributions 
of regularly varying tails and (sub)exponential tails. However, we are not basically concerned with 
any distribution with a superexponential tail, e.g., a multivariate normal distribution. The details 
of the FCLT in that case remain for a future study. 

The remainder of the paper is organized as follows. First, in Section 2 we provide a formal 
definition of the subgraph counting process. Section 3 gives an overview of what was shown in the 
previous work [20| and what will be shown in this paper. Subsequently, in Section 4 we focus on 
the case in which the underlying density has a regularly varying tail, including power-law tails, and 
prove the required FCLT. We also investigate the properties of the limiting Gaussian processes, 
in particular, in terms of self-similarity and sample path continuity. In Section 5, we do the same 
when the underlying density has an exponentially decaying tail. To distinguish densities via their 
tail properties, we need basic tools in EVT. In essence, the properties of the limiting Gaussian 
processes are determined by how rapidly Rn grows to infinity, as well as how rapidly the tail of / 
decays. Finally, Section 6 carefully examines both cores and weak cores for a large class of densities. 

Before commencing the main body of the paper, we remark that all the random points in this 
paper are assumed to be generated by an inhomogeneous Poisson point process on with intensity 
nf. In our opinion, the FCLT in the main theorem can be carried over to a usual i.i.d. random 
sample setup by a standard “de-Poissonization” argument; see Section 2.5 in |21] . This is, however, 
a little more technical and challenging, and therefore, we decided to concentrate on the simpler setup 
of an inhomogeneous Poisson point process. Furthermore we consider only spherically symmetric 
distributions. Although the spherical symmetry assumption is far from being crucial, we adopt it 
to avoid unnecessary technicalities. 


2. Subgraph Counting Process 

Let (Xj, z > I) be i.i.d. M'^-valued random variables with spherically symmetric probability 
density /. Given a Poisson random variable Nn with mean n, independent of (Xj, z > 1), denote 
by Rn = {Xi,X 2 ,... ,X 7 v„} a Poisson point process with \Vn\ ■= X„. We choose a positive integer 
k, which remains fixed hereafter. We take k > 2, unless otherwise stated, because many of the 
functions and objects to follow are degenerate in the case of A: = 1. 

Let F be a fixed connected graph of k vertices and G represent a geometric graph; = denotes 
graph isomorphism. We define 

/z(xi, ...,Xk):= l{G({xi,... ,Xfc}, I) ^ F} , xi,..., Xfc € M"*. 

Next, we define a collection of indicators (hj, t > 0) by 

(2.1) ht{xi,.. .,Xk) := h{xi/t,.. .,Xk/t) = 1{G({xi, ... ,Xfc}, t) = r} , 
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from which one can capture the manner in which a geometric graph dynamicaiiy evoives as the 
threshoid radius t varies. Note, in particuiar, that , Xk) = h{xi ,..., Xk)- 

Cieariy ht is shift invariant: 

(2.2) ht{xi,...,xk) = ht{xi + y,...,xk + y), xi,..., y G M"*, 
and, further, 

(2.3) /it(0, xi,..., Xfc-i) = 0 if ||xi|| >/ct for some i = 1,..., A: — 1. 

The iatter condition impiies that /it(xi,... ,Xfc) = 1 oniy when aii the points xi,... ,Xfc are ciose 
enough to each other. 

Moreover ht can be decomposed as foiiows. Suppose that T has k vertices and j edges for some 
j G {/c — 1,...,A;(A: — l)/2}. Letting be a set of connected graphs of k vertices and I edges (up 
to graph isomorphism), define for xi,... ,Xfc G 

k(k-l)/2 

hf{xi,...,xk)--=ht{xi,...,xk)+ ^ l{G({xi,... ,Xfc},t) = r'l , 

i=j+l 

k(k-l)/2 

V(xi,... ,Xfc) := ^ ^ l{G({xi,...,Xfc},t) =r'} . 

t=j+l T'&Ai 

Note that /i)^(xi,... ,Xfc) = 1 if and oniy if a geometric graph G({xi,... ,Xk},t) either coincides 
with r (up to graph isomorphism) or has more than j edges, whiie h^{xi,... ,Xk) = 1 oniy when 
G({xi, ..., Xfc}, t) has more than j edges. It is then eiementary to check that hf are both indicators, 
taking vaiues 0 or 1, and satisfying, for aii xi,..., x^ G and 0 < s < t. 


(2.4) 

ht{xi, ...,Xk) = hj+(xi,... ,Xfc) 

- /l7(xi,...,Xfc), 


(2.5) 

/l+(xi,...,Xfc) < /i 4 +(xi,.... 

Xk) , 



h~{xi,...,Xk) < V(xi,..., 

Xk) ■ 


(2.6) 

hf{0,xi,.. . ,Xfc_i) = 0 if II 

Xi\\ > kt for some i = 1,.. 

. ,k — 1 


In addition, since ht is an indicator, it is aiways the case that 

ht{xi,...,Xk) < hf{xi,...,Xk). 

The objective of this study is to estabiish a functionai centrai iimit theorem (FCLT) of the 
subgraph counting process defined by 

(2.7) Gn{t)-.= Y, ^t{y)l{m{y)>Rn], t>0, 

ycVr. 

where ht is given in (12.11) . m(xi,... ,Xk) = mini<j<fc ||xi||, Xj G and {Rn, n > 1) is a property 
chosen normaiizing sequence. Note that ()2.7p counts the number of subgraphs in G{'Pn,t) isomor¬ 
phic to r that iie compieteiy outside of B{0,Rn). More concrete definitions of (Rn) are given in 
the subsequent sections, where the sequence is shown to be dependent on the tail decay rate of /. 

3. Annuli Structure 

The objective of this short section is to clarify what is already known and what is new in 
this paper. Without any real loss of generality, we will do this via two simple examples, one of 
which treats a power-law density and the other a density with a (sub)exponential tail. Before this, 
however, we introduce two important notions. 
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Definition 3.1. m) Given an inhomogeneous Poisson point process Vn in with a spherically 
symmetric density /, a centered ball B{0, Rn), with Rn —)• oo, is called a core if 

(3.1) B{0,Rn)C U B{X,1). 

XG'PrLr]B[O^Rn) 

In other words, a core is a centered ball in which random points are densely scattered, so that 
placing unit balls around them covers the ball itself. We usually wish to seek the largest possible 
value of Rn such that occurs asymptotically with probability 1. A related notion, the weak 
core, plays a more decisive role in characterizing the FCLT proven in this paper. It is shown later 
that a weak core is generally larger but close in size to a core of maximum size. 

Definition 3.2. Let / be a spherically symmetric density on and ei = (1,0,... ,0) € A 
weak core is a centered ball B{0, R^n'^) such that nf{R^'^ei) —>■ 1 as re —>■ oo. 

Example 3.3. Consider the power-law density 

(3.2) fix) = C/{l + \\x\\^) , 

for some a > d and normalizing constant C. Using this density, we see how random geometric 
graphs are formed in all of M'^. First, according to [T], there exists a sequence Rn'^ ~ constant x 
(re/logre)^/“, re —)■ oo such that, if Rn < Rn\ (j3.ip occurs asymptotically with probability 1. In 
addition, as for the radius of a weak core, it suffices to take rI^'^ = (Cre)^/". Although rI^^ grows 
faster than Rn , they are seen to be “close” to each other in the sense that they have the same 
regular variation exponent, 1/a. 

Beyond a weak core, however, the formation of random geometric graphs drastically varies. In 
fact, the exterior of a weak core can be divided into annuli of different radii, at which many isolated 
subgraphs of finite vertices are asymptotically placed in a specific fashion. To be more precise, let 
us fix connected graphs F^ with k vertices for k = 2,3,... and let 

which in turn implies that R^^'^ ^ • <C R^}^ -C R^^li „ ^ ^ s-iid 

^ 1 > n ^ oo . 

Under this circumstance, [20] considered the subgraph counts given by 

(3.3) ^ l{G(T,t) = F4xl{rre(T)><i}, 

ycVn 

and showed that (j3.3p weakly converges to a Poisson distribution for each fixed t. To be more 
specific on the geometric side, let Ann(A', L) be an annulus with inner radius K and outer radius 
L. Then, we have, in an asymptotic sense, 

• Outside B{^,R^^^, there are finitely many graphs isomorphic to F 2 , but none isomorphic 
to F 3 ,F 4 , .... 

• Outside i?( 0 ,i? 3 ^^), equivalently inside Ann(i? 3 ^^, , there are infinitely many graphs 

isomorphic to F 2 and finitely many graphs isomorphic to F 3 , but none isomorphic to 

F4,F5, .... 

In general. 
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Figure 1. Layered structure of random geometric graphs. For the density (13.21) . Rn^ and 
are regularly varying sequences with exponent a~^. R/^l^ is also a regularly varying sequence with 
exponent (a—d/k)~^. We study the FCLT for (|3.4I) in three different regimes, i.e., (i) nf{R„ei) —>■ 0, 
(ii) n/(i?„ei) —>■ ^ € (0,oo), and (in) n/(i?„ei) —>■ oo. In relation to other radii, they are 
respectively equivalent to (i) r!'™^ R„ ~ and (Hi) R„ -C 


• Outside B (O, , equivalently inside Ann(i?^^^, „), there are infinitely many graphs 

isomorphic to r 2 ,..., F^-i and finitely many graphs isomorphic to F^, but none isomorphic 
to Ffc+i, rfc+ 2 , • • • etc. 

Section [5] of the current paper considers the subgraph counts of the form 

(3.4) ^ ^ Ffc} X l{m(3^) > Rn] , 

ycVr. 

where (Rn) satisfies 

(3.5) R!^f (Rnei)^ ^ oo , n —>• oo , 

in which case, Rn <C R^^n- As a consequence of (j3.5p . we may naturally anticipate that a FCLT 
governs the asymptotic behavior of (j3.4p . Since (Rn) satisfying (j.S.5p shows a slower divergence 
rate than (R^j^lj, i.e., Rn/R^^\ 0, we may expect that infinitely many subgraphs isomorphic to 

Ffc appear asymptotically outside B({),Rn). This in turn implies that, instead of a Poisson limit 
theorem, the FCLT governs the limiting behavior of the subgraph counting process. 

As the analog of the setup for (II.ip . when deriving an FCLT, the behavior of (Id.dh splits into 
three different regimes: 

(i) nf(Rnei) ^ , (ii) n/(i?nei) ^ ^ G (0, oo), (in) nf(Rnei) ^ oo . 

Specihcally, if nf(Rnei) —0 (i.e., B(0,Rn) contains a weak core), many isolated components of 
subgraphs isomorphic to T^ are distributed outside B(0,Rn). If nf(Rnei) -H- ^ £ (0,oo) (i.e., 
B(0,Rn) agrees with a weak core), the subgraphs isomorphic to T^ outside B(0,Rn) begin to be 
connected to one another. In particular, observing that limfc_).oo R^kn ~ 

• Outside of i?(0, there are infinitely many graphs isomorphic to T^ for every j = 

2,3,.... 

If nf(Rnei) —>• oo (i.e., B(0,Rn) is contained in a weak core), the subgraphs isomorphic to T^ 
outside B(0,Rn) are further increasingly connected and form a large component. 

In Section 01 we will see that the nature of the FCLT, including the normalizing constants and 
the properties of the limiting Gaussian processes, differs according to which regime one considers. 
Combing the results on the FCLT and the Poissonian results in [20] , we obtain a complete picture 
of the annuli structure formed by heavy tailed random variables. 
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Example 3.4. Next, we turn to a density with a (sub)exponential tail 

fix) = X G 0 < r < 1. 

for which the radius of a maximum core is given by 

= (r log n — T log log(T log + constant) 

see [1] and [20]. Obviously, one can take = (r log n + r log C) As in the previous example, 
the exterior of a weak core is characterized by the same kind of layer structure, for which the 

(p) 

description in Figure 1 applies, except for the change in the values of Letting 

= (r logn + - r) log(Tlogn) + r log C) , 

we have, in an asymptotic sense, <C • • • <?; R^^ ^ n ^ 

n^oo. 

Then, it was shown in [20| that (13.31) converges weakly to a Poisson distribution for each fixed t. 

In Section [5] of this paper, taking [Rn) such that vf R'^'^f {Rnei)^ —>■ oo, we establish a FCLT 
for the subgraph counting process (|2.7I) . To this end, our argument has to be split, once again, into 
the three different regimes: 

(f) nfiRnei) 0 , {ii) nfiRnCi) ^ G (0, oo), iiii) nfiRnCi) oo . 

As in the last example, three different Gaussian limits may appear depending on the regime. 
This completes the full description of the annuli structure formed by random variables with an 
exponentially decaying tail, when combined with the Poisson limit theorems in m- 


4. Heavy Tail Case 


4.1. The Setup. In this section, we explore the case in which the underlying density / on has a 
heavy tail under a more general setup than that in Example 13.31 Let Sd-i be a {d— l)-dimensional 
unit sphere in M'^. We assume that the density has a regularly varying tail (at infinity) in the sense 
that for any 6 G Sd-i (equivalently, for some 9 G Sd-i because of the spherical symmetry of /), 
and for some a > d, 


lim 

r—¥oo 


fjrte) 

fire) 


= t-“ 


for every f > 0 . 


Denoting by RV-a a collection of regularly varying functions (at infinity) of exponent —a, the 
above is written as 


(4.1) / G RV-a ■ 

Clearly, a power-law density in Example 13.31 satisfies (|4.1I) . Let A: > 2 be an integer that remains 
fixed throughout this section. We remark that many of the functions and objects are dependent 
on k, but the dependence may not be stipulated by subscripts (or superscripts). Choosing the 
sequence Rn oo so that 

(4.2) R!^f iRnCi)^ oo as n —oo , 

we consider the subgraph counting process given in (|2.7I) . whose behavior is, as argued in Example 
13.31 expected to be governed by a FCLT. 
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The scaling constants for the FCLT, denoted by r„, are shown to depend on the limit valne of 
nf{Rnei) as n —>■ oo. More precisely, we take 

'n>^RifiRnei)>^ if n/(i?„ei) ^ 0 , 

Ri if n/(i?„ei)c G (0,oo), 

n^'^-^R^f{Rneif^-^ if n/(i?„ei) oo . 


(4.3) 


Tn ■ = 


The reason for which we need three different normalizations is deeply related to the connectivity 
of a random geometric graph. To explain this, we need the notion of a weak core] see Definition 
[Q for the formal definition. The main point is that the density of random points between the 
ontside and inside of a weak core is completely different. In essence, random points inside a weak 
core are highly densely scattered, and the corresponding random geometric graph forms a single 
giant component. Beyond a weak core, however, random points are distributed less densely, and 
as a result, we observe many isolated geometric graphs of smaller size. This disparity between the 
outside and inside of a weak core requires different normalizations in (Tn)- In Section [6l a more 
detailed study in this direction is presented. 


4.2. Limiting Gaussian Processes and the FCLT. We introduce a family of Gaussian pro¬ 
cesses which function as the building blocks for the limiting Gaussian processes in the FGLT. For 
£ = 1, . . . , k, let 

^ l\{{k - {a{2k - tj - d) ' 

where Sd_i is a surface area of the {d — l)-dimensional unit sphere in W^. 

For £ = 2,... ,k, write Xg for the Lebesgue measure on and denote by Gg a Gaussian 

BiXi-noise, such that 

G,iA)r^M{0,BAi{A)) 

for measurable sets A C with X(,{A) < oo, and if ^ n B = 0, then G(,{A) and Gi{B) are 

independent. For ^ = 1, we define Gi as a Gaussian random variable with zero mean and variance 
Bi. We assume that Gi,... ,Gk are independent. 

For £ = 2,..., k — 1, we define Gaussian processes t > O) by 


In addition, if £ 


and if ^ = 1, set 


Viit) ■= f [ ht{0,y,z)dzGi{dy), t>0. 
k, define 

Vk{t) := [ ht{0,y)Gk{dy), 

J(IRd)fc-l 


Viit) ■= [ ht{0,z) dzGi = [ h{0,z)dzGi. 

Note that Vi is a degenerate Gaussian process with deterministic sample paths. These processes 
later turn out to be the building blocks of the weak limits in the main theorem. 

The covariance function of the process Vf is given by 


(4.4) 


Li{t,s) -ElViimis)} 


Bi 




dz2 

,k-i 



dzi ht( 0 ,y,zi) /i^(0,y,Z2), 

,k-i 


(if £ = k, we take Zj = 0, f = 1, 2, and if ^ = 1, we set y = 0). 


f, s > 0 
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Using the decomposition (j2.4|) . we can express as the difference between two Gaussian pro¬ 
cesses; that is, for ^ = 2,..., A: — 1, 

Ve{t)=[ [ ht{0,y,z)dzGe{dy)- [ [ {0,y,z) dzGi{dy) 

:=V+it)-Vf{t). 

The same decomposition is feasible in an analogous manner for Vi and V^. 

The following proposition shows that the processes and can be represented as a time- 
changed Brownian motion. 

Proposition 4.1. The process can be expressed as 

(Ufc+(i), i > 0 ) = t > o), 

where B is the standard Brownian motion, and := B^ h'^{0,y)dy. 

Replacing with := B^ f^j^djk-i h~{0,y)dy, we obtain the same statement for V^. 

Proof. It is enough to verify that the covariance functions on both sides coincide. It follows from 
(j2.5D that for 0 < s < t, 

[ ht{0,y)ht{0,y)dy 
= K{B{K+ t^^'‘-^^)B{K+ 

□ 

We also claim that the process is self-similar and has a.s. Holder continuous sample paths. 
Recall that a stochastic process [X{t), t > O) is said to be self-similar with exponent H if 

{X{cti), i = l,...,k) = {c^X{ti), i = l,...,k) 
for any c > 0, ti,..., > 0, and fc > 1. 

Proposition 4.2. (i) For I = 1,..., A:, the process 'Vi is self similar with exponent H = d(2k — 
£-l)/2. 

(ii) For £ = 1,... ,k and every T > 0, [Vi{t), 0 < t < T) has a modification, the sample paths of 
which are Holder continuous of any order in [0,1/2). 

Proof. We can immediately prove (z) by the scaling property 

Li{ct, cs) = Li{t, s), t,s > 0 , c > 0 . 

As for (n), the statement is obvious for £ = 1 or £ = A;; therefore, we take £ G {2,..., A: — 1}. By 
Gaussianity, 

m 

(4.5) E{(U,(t)-'^/(s))'™} m = l,2,... 

i=l 

We now show that there exists a constant C > 0, which depends on T, such that 

(4.6) E|(W(t)-W(s))^} <C(t-s) for all 0 < s < t < T. 
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By virtue of the decomposition , showing (j4.6|) for each of and suffices. We 

handle only, since can be treated in the same manner. We have 

E|(y/(t) - / dy [ dz 2 [ dzi{/i^(0, y, zi) -/i+(0, y, zi)} 

J JiRdy-i J(Rd)k-e 

X {hi (0, y, Z 2 ) - hi (0, y, Z 2 )} . 

Because of (12.6p . the above integral is not altered if the integral domain is restricted to x 

^ (B(0,feT))^ In addition, by (12.51) . there exist constants Ci,C 2 > 0, both depending 
on T, such that 

M (t) - {s)y] < Cl [ [ [hi{ 0 ,y,z) - hi{ 0 ,y,z)}dzdy 

J JiRdy-i J(Rdy-i 

= Ci [ [ /i+(0, y, z)dzdy 

J{Rdy-t 

<C 2 / ^"*"(0, y, z)dzdy {t — s) for all 0 < s < t < T, 

J{Rdy-t 

which verihes (|4.6p . 

Combining (|4.5p and (|4.6I) . we have that for some C 3 > 0, 

< Csit - sr for all 0 < s < t < T. 

It now follows from the Kolmogorov continuity theorem that there exists a modihcation of (V^(t), 0 < 
t the sample paths of which are Holder continuous of any order in [O, {m — l)/(2m)). Since 

m is arbitrary, we are done by letting m ^ 00 . □ 

We are now ready to state the FCLT for the subgraph counting process, suitably scaled and 
centered in such a way that 

Xn{t) = T-C^{Gn{t) - E{G„(t)}) , t > 0 . 

In the following, ^ denotes weak convergence. All weak convergence hereafter are in the space 
X>[0, cx)) of right-continuous functions with left limits. The proof of the theorem is deferred to 
Section 17.11 

Theorem 4.3. (i) If nf{Rnei) ^0 as n ^ 00 , then 

[Xn{t), t > 0 ) => (Vk{t), t > 0 ) in T>[0, 00 ). 

(ii) If nf{Rnei) —^ G (0, 00 ) as n ^ 00 , then 

{Xn{t), t > 0 ) ^ (j 2 f’'~^Vi{t), t > 0^ in P[0,oo). 

(Hi) If nf{Rnei) —>■ 00 as n ^ 00 , then 

[Xn{t), t > 0 ) ^ (^ 1(^)5 i > 0 ) in ^[O, 00 ). 

The processes Vi,..., can be viewed as the building blocks of the limiting Gaussian processes; 
however, how many and which ones contribute to the limit depends on whether the ball B{ 0 ,Rn) 
covers a weak core or not. If B{0,Rn) covers a weak core, equivalently, nf{Rnei) 0, then 

is the only process remaining in the limit. Although, as seen in Proposition 14.11 is generally 

represented as the difference in two time-changed Brownian motions, it can be denoted as a single 
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time-changed Brownian motion when ht is increasing in t, i.e., hs{y) < ht{y) for all 0 < s < t, 3^ € 
This is the case when T is a complete graph, in which case the negative part hj' is identically 
zero. In contrast, the process Vi, a degenerate Gaussian process with deterministic sample paths, 
only appears in the limit when B{0, Rn) is contained in a weak core, i.e., nf{Rnei) — )• oo. Finally, 
if 5(0, Rn) agrees with a weak core (up to multiplicative constants), all of the processes Vi,..., 
contribute to the limit. Interestingly, only in this case, do the weak limits become non-self-similar. 


5. Exponentially Decaying Tail Case 


5.1. The Setup. This section develops the FCLT of the subgraph counting process suitably scaled 
and centered, when the underlying density on possesses an exponentially decaying tail. Typically, 
in the spirit of extreme value theory, a class of multivariate densities with exponentially decaying 
tails can be formulated by the so-called von Mises functions. See for example, [3] and [1]. In 
particular, in the one-dimensional case (d = 1), the von Mises function plays a decisive role in the 
characterization of the max-domain of attraction of the Gumbel law. See Proposition 1.4 in |22j . 
We assume that the density / on is given by 

(5.1) /(x) = L(||x||) exp{-'0(||a:||)} , x G 

Here, ^|J : M+ —)• M is a function of C^-class and is referred to as a von Mises function, so that 

(5.2) > 0, ip{z) oo, (l/ip'y{z) 0 


as z ^ Zoo S (0, oo]. In this paper, we restrict ourselves to an unbounded support of the density, 
i.e., Zoo = OO- For notational ease, we introduce the function a(z) = l/'0'(z), z > 0. Since a'{z) —)• 0 
as z ^ oo, the Cesaro mean of a' converges as well: 

/ ^ aiz) 1 r ,, , , 

(5.3) - = - a {r)dr —>■ 0 , as z —)■ oo . 

^ ^ Jo 

Suppose that a measurable function L : M_|_ ^ M_|_ is flat for a, that is. 


(5.4) 


L(t -|- a{t)v) ^ 

W) 


as t —>■ oo uniformly on bounded u-sets. 


This condition implies that L behaves as a constant locally in the tail of /, and thus, only if plays 
a dominant role in the characterization of the tail of /. Here, we need to put an extra technical 
condition on L. Namely, there exist 7 > 0, zq > 0; and C >1 such that 


(5.5) 


L{zt) 

L{z) 


< CB 


for a\\ t > 1, z > zq . 


Since L is negligible in the tail of /, it seems reasonable to classify the density (|5.ip in terms of the 
limit of a. If a{z) oo as z ^ oo, we say that / belongs to a class of densities with subexponential 
tail, because the tail of / decays more slowly than that of an exponential distribution. Conversely, if 
a{z) ^ 0 as 2 ; —>■ 00 , / is said to have a superexponential tail, and if a{z) ^ c G (0, 00 ), we say that 
/ has an exponential tail. To be more specific about the difference in tail behaviors, let us consider a 
slightly more general example than that in Example 13.41 for which /(x) = L(||x||) exp| — ||x||'^/t}, 
r > 0, X G M'^. Clearly, the parameter r is associated with the speed at which / vanishes in the 
tail. Observe that a{z) = z^~'^ —^-ooasa;—>'Ooif0<T<l, and therefore in this case, / has a 
subexponential tail. If r > 1, a{z) decreases to 0, in which case / has a superexponential tail. 

An important assumption throughout most of this study is that there exists c G (0, 00 ] such that 


( 5 . 6 ) 


a(z) —>■ c as 2; —>■ 00. 
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In view of the classification described above, (|5.6jl eliminates the possibility of densities with su¬ 
perexponential tail. As discovered in [ 20 ] and [T], random points drawn from a superexponential 
law hardly form isolated geometric graphs outside a core, whereas random points coming from a 
subexponential law do constitute a layer of isolated geometric graphs outside a core. Accordingly, 
it is highly likely that the nature of the FCLT differs according to whether the underlying density 
has a super exponential or a sub exponential tail. The present work focuses on the (sub)exponential 
tail case, and more detailed studies on a superexponential tail case remain for future work. 

To realize a more formal set up, let A; > 2 be an integer, which remains fixed for the remainder 
of this section; however, once again, note that many of the functions and objects are implicitly 
dependent on k. Dehne the sequence Rn —>■ oo, so that 

(5.7) n^a{Rn)Rn~^ f{Rnei)^ ^ oo , n ^ oo . 

(p) 

Defining an alternative sequence i?). „ —>■ oo for which 

^ 1 , n ^ 

(P) 

the subgraph counting process using R^^!^ is known to weakly converge to a Poisson distribution; 

see [20]. Since Rn in (j5.7p grows more slowly than R^]^, i.e., Rn/R^^\ —>■ 0, we may expect that an 
FCLT plays a decisive role in the asymptotic behavior of a subgraph counting process. 

As in the last section, we now want to recall the notion of a weak core. Let R^h'^ ^ oo be a 
sequence such that nf{Rl^^ei) —>■ 1 as n —>■ oo. Then, we say that a ball B{0,rI^^) is a weak 
core. We have to change, once again, the scaling constants of the FCLT, depending on whether 
B{0, Rn) covers a weak core or not. More specifically, we dehne 


(5.8) 


'^n • — 


'n’^a{Rn)Rt^fiRnei)^ 

< a{Rn)Ri~^ 

n^>^-^aiRn)Rt^fiRnei)^^-^ 


if nf{Rnei) 0 , 

if nfiRnCi) C € ( 0 ,oo), 

if nfiRnCi) oo . 


5.2. Limiting Gaussian Processes and the FCLT. The objective of this subsection is to for¬ 
mulate the limiting Gaussian processes and the FCLT. Let 

(5.9) Di = - — -^ = 1,...,A:, 

and let Hi be a Gaussian /r^-noise, where the m ioi £ = 2,..., k, satisfy 
^^{dpdY) = Ei=i(ei.yi> 

X l{p + c~^{ei,yi) >0, z = 1,... - 1 } dpdy, p > 0, y € (M'^)^"\ 

and 

pi{dp) = Di e~^dp, p > 0. 

Assume that Hi, ..., are independent. 

We now dehne a collection of Gaussian processes needed for the construction of the limits in the 
FCLT. For i = 2,... ,k, we dehne 

We{t):=[ [ 

X l{ p + c~^ {ei, Zi) > 0, i = 1,... ,k - £} ht(0,y,z) dzH^(dpdy), 
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and, accordingly, 

Wi{t) := [ [ 

X l| /? + c~^{ei,Zi) >0, i = 1,... ,k — ht{0,z) dz Hi (dp), 

Wk{t) := f ht{0,y) Hkidpdy). 

J[0,oo)x(Rd)'=-i 

As we did in Section [4.21 by the decomposition ht = hf — h~[, one can write the process as the 
corresponding difference for ^ = 1 ,..., A:. 

It is easy to compute the covariance function of W^. We have, for i = 1,... ,k and t,s >0, 


(5.10) Mi{t,s):=E{We{t)Weis)} 


= D, 


roo r 

Jo J(M 


lo J{Rd-) 2 k-t-i 

1/„ „.\^n „• _ 1 itu 0 l,W I 


xl{p + c ^(ei,yi)>0, i = 1,... ,2/i:l}/r)j(0,y)(iy d/ 9 , 

where 

(5-11) h[^l ( 0 , 2 / 1 ,..., y 2 k-l-i) ■■= ht{0, 2 / 1 ,..., 2 /fc-i) ^s(0, 2 /i, • • •, 2 /£-i, 2 /fc, ■ ■ ■, 2 / 2 fc-£-i), 
and, in particular, we set 


^s(0, 2 / 1 , • • •, 2/£-i, 2/fc, ■ ■ ■ , y 2 k-e-i) 


hs{0,yk,---,y2k-2) if^ = l, 

/is(0,2/1,... ,2/fc-i) if^ = A:. 


It is important to note that if a{z) —>■ oo as z —>■ oo, then Mi coincides with Li given in (14.411 up 
to multiplicative factors, i.e., 

M£{t, s) = [a — d{2k — s), t,s > 0. 


This in turn implies that 

W£={a-d{2k-i)-^y^Wi, 

in which case, there is nothing to explore here, because the properties of have already been 
studied in Section 14.21 

In contrast, if a{z) —c G (0, oo) as z —>■ oo, then Mi does not directly relate to Li as above, and, 
consequently, the process exhibits properties different to those of V^. For example, although 
one may anticipate, as the analog of the process Vi, that Wi is a degenerate Gaussian process, 
this is no longer the case. 


Proposition 5.1. Suppose that a{z) — )• c G (0, oo) as z ^ oo. 

(i) Wi is a non-degenerate Gaussian process. 

{ii) For 1= 1,..., A;, is non-self-similar. 

Proof. If a{z) —>■ c G (0, oo) as z —>■ oo, then Mi{t,s) cannot be decomposed into a function of t 
and a function of s, and therefore, Wi is non-degenerate. 

As for (li). Ml does not match Li at all and it loses the scale invariance, meaning that is 
non-self-similar. □ 


Similarly to Proposition 14.11 however, the process W/j(= — W^) can be denoted in law as 

the difference between two time-changed Brownian motions, regardless of whether a{z) —)• oo or 
a{z) ^ c G (0, oo) as oo. Furthermore, the sample paths of are Holder continuous. 
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Proposition 5.2. Irrespective of the limit of a, the following two results hold. 

{i) The process can be represented in law as 

(M 4 +(t), t>0) = (nf [ h^{0,y) Hkidpdy)], t > o) , 

where B is the standard Brownian motion. 

The same statement holds for W^, by replacing hi^ with h'f. 

(ii) For i = 1,... ,k, and every T > 0, {Wi{t), 0 < t < T) has a modification, the sample paths of 
which are Holder continuous of any order in [0,1/2). 

Proof. The proof of (i) is very similar to that in Proposition 14.11 so we omit it. The proof of (ii) 
is analogous to that in Proposition 14.21 (ii)] we have only to show that for some C > 0, 

^[{Wi{t) -Wi{s)f"^ <C{t-s) for all 0 < s < t < T. 

Because of the decomposition — W^, it suffices to prove the above for each and 

W^. We check only the case of W^. We see that 



X (h)^(0 ,y,z) - /i+(0 ,y,z)) dzj pi{dpdy) 
<D,Bfi^¥.[{V+{t)-V+{.s)f]. 

The rest of the argument is completely the same as Proposition 14.21 [ii). □ 

Now, we can state the FCLT of the centered and scaled subgraph counting process 

Xn{t) = T-^/^{Gn{t) -nGn{t)}) , t >0, 

where the normalizing sequence [Rn) satisfies (15.71) and (t„) is defined in (15.81) . Interestingly, if 
/ has a subexponential tail, i.e., a[z) —>■ oo, then the limiting Gaussian processes in the theorem 
below completely coincide (up to multiplicative constants) with those in Theorem 14.31 When / 
has an exponential tail, i.e., a[z) ^ c € (0, oo), the limiting Gaussian processes are essentially 
different from those in Theorem 14.31 The proof of the theorem is presented in Section 17.21 For 
the reader’s convenience, we summarize in Tables 1 and 2 the properties of the limiting Gaussian 
processes in Theorems 14.31 and 15.31 These tables indicate that the limiting Gaussian processes are 
somewhat special when / has an exponential tail. For example, in this case, the limits always lose 
self-similarity, regardless of the asymptotics of nf[Rnei), whereas, in the regularly varying or the 
subexponential tail case, the self-similarity is lost only when nf[Rnei) converges to a positive and 
finite constant. Furthermore, when nf[Rnei) oo, a non-degenerate limit appears only in the 
exponential tail case. 

Theorem 5.3. Assume that the density ()5.ip satisfies (15.2p . ()5.4I) . (15.51) . and ()5.6p . 

[i) If nf[Rnei) 0 as n ^ oo, then 

[Xn{t), t > 0) ^ {Wk[t), t > 0) in T)f), oo). 
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(m) If nf{Rnei) f, € (0, oo) as n ^ oo, then 


{Xn{t), t>0)^ t>0^ in V[0,oo) 


(Hi) If nf{Rnei) —>■ oo as n ^ oo, then 

[Xn{t), t > O) (Wi{t), i > O) in "Dp, oo). 



nf{Rnei) 0 

nf{Rnei) ^ 

n/(R„ei) oo 

Regularly varying tail 

d{k - l)/2 

Non-SS 

d{k — 1) 

Subexponential tail 

d{k-l)/2 

Non-SS 

d{k — 1) 

Exponential tail 

Non-SS 

Non-SS 

Non-SS 


Table 1. Self-similarity exponents of the limiting Gaussian processes. Non-SS 
means that the process is non-self-similar. A zero limit of nf{Rnei) is equivalent 
to the case in which a ball B{0,Rn) contains a weak core, and nf{Rnei) —)■ oo if 
and only if B{0,Rn) is contained in a weak core. If nf{Rnei) G (0,oo), then 
B{0,Rn) agrees with a weak core (up to multiplicative constants). 



nf{Rnei) 0 

nf{Rnei) ^ 

nf{Rnei) —>■ oo 

Regularly varying tail 

Difference of time-changed 
Brownian motions 

New 

Degenerate Gaussian 
process 

Subexponential tail 

Difference of time-changed 
Brownian motions 

New 

Degenerate Gaussian 
process 

Exponential tail 

Difference of time-changed 
Brownian motions 

New 

New 


Table 2. Representation results on the limiting Gaussian processes. “New” implies 
that the limit constitutes a new class of Gaussian processes. 


6. Graph Connectivity in Weak Core 

We start this section by recalling the weak core, which was defined as a centered ball i?(0, 
such that ^ 1 as n ^ oo. In addition, we need the relevant notion, the core, which was 
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defined in Definition 13.11 Recall that, given a Poisson point process Vn on a core is a centered 
ball B{0,Rn) such that 

(6.1) i?(0,R„)C (j 

X G7^nnS(0,i?ri) 

In the following, we seek the largest possible sequence Rn —>■ oo such that the event ()6.ip occurs 
asymptotically with probability 1, and subsequently, it is shown that the largest possible core and 
a weak core are “close” in size. However, the degree of this closeness depends on the tail of an 
underlying density /, and therefore, we divide the argument into two cases. 

We first assume that the density / on is spherically symmetric and has a regularly varying 
tail, as in (14.ip . For increased clarity, we place an extra condition that p{r) := f{rei) is eventually 
non-increasing in r, that is, p is non-increasing on (ro, oo) for some large vq > 0. In this case, the 
radius of a weak core is, clearly, given by 

(6.2) R^^^ = (n) := infjs ; (s) > nj. 

Proposition 6.1. Suppose that p € RV-a for some a > d and p is eventually non-inereasing. 
Define 

/ dm 

p J \ log n — 82 log log n 

with G (0, Q;/(2'^(i'^/^+^)) and 82 G (0,1). If Rn < Rn\ then 



(6.4) P| B{0,Rn) C B{X,1) 1^1, n —>■ 00. 

Furthermore, the sequences {Rn^) in (I6.3h and {R^^) in (16.2h are both regularly varying sequences 
with exponent 1/a, and 


(6.5) 


p(c) 

rCn 


R 


(w) 


81 


1/0 


0 , n —>■ 00 . 


^ log n — 82 log log n ^ 

One can obtain a parallel result when the underlying density has an exponentially decaying tail, 
as in m- We simplify the situation a bit by assuming 


(6.6) f{x) = Cexp{-fi{\\x\\)}, xGM'^, 

where C is a normalizing constant and : M+ ^ M is of O^-class and satisfies G RV^ (at infinity) 
for some u > 0 and > 0. It should be noted that we are permitting the case u > 1, implying that, 
unlike in the previous section, we do not rule out densities with superexponential tail. Evidently, 
the radius of a weak core is given by 

(6.7) ='i/^(logn -I- logC). 

Proposition 6.2. Assume that a probability density f on is given by (16.6p . Define 

(6.8) (log n — log log log n — (5i — ^2), 

where (5i = dlog 2 — log u -|- (1 -|- d/ 2 ) log d — log C and 82 > 0. If Rn < Rn \ then 


B{0,Rn)C U B{X,1) 

^Hn ) 


(6.9) 


P 


n —)■ 00 . 
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Furthermore, the sequences {Rn^) in (|6.8I) and in (j6.7p are close in size in the sense of 

log log log re + (^1 + (52 + log C\ 
log Cn ) 


( 6 . 10 ) 


p(c) 


- 1 - 


0 , 


oo. 


The following result is needed as preparation for the proof of these propositions. The proof 
may be obtained by slightly modifying the proof of Theorem 2.1 in [T], but so that this paper is 
self-contained, we repeat the argument. 


Lemma 6.3. Given a spherically symmetric density f ore suppose that p{r) = f{rei) is eventu¬ 
ally non-increasing. Let g = l/(2(i^/^). Suppose, in addition, that there exists a sequence Rn oo 
such that dlogRn — g^nf{Rnei) ——(X) as n ^ oo. Then, 


(6.11) Pj B{0,Rn)C B{X,1) 1 —)• 1, re ^ oo. 

\ X&Vr.r\B{0,Rn) ) 

Proof. For p >0, let Q{p) be a collection of cubes with grid g that are contained in B{t), p). Then, 
{QnK/0 for alige Q(p)} C {b(0,p) C IJ B{x,l)j 

xeVnr\B{o,p) 

for all /o > 0 and re > 1. It now suffices to show that 

P( g n T’n = 0 for some Q G Q{Rn) ) ^ 0 , re ^ oo . 

This probability is estimated from above by 

^ P(gnPn = 0 )= ^ exp|-re / /(x)dx| 


Q^Q{Rn) 


Q^Q{Rn) 


< ^ exp{-re//(i2nei)} < 5 '^^nexp{-/re/(i?„ei)} . 

QeQ{Rn) 

At the first inequality, we used the fact that p is eventually non-increasing. Clearly, the rightmost 
term vanishes as re ^ oo. □ 


Proof, {proof of Provosition [h!7]l Observe that the assumption p G RV-a implies (1/p)^ G RVii^, 
e.g.. Proposition 2.6 (v) in [23]. Thus, (|6.5I) readily follows from the uniform convergence of 
regularly varying functions; see Proposition 2.4 in [23] . By Lemma 16.31 it suffices to verify that 
dlogiZn^ - g‘^nf{Rn^ei) —>■ —oo as re ^ oo. Since 0 < (52 < 1, we have 


dlogi?(f) < d 


log ( - 
P 


6 in 


log re — (52 log log re 
X (log (5i -I- log re - 62 log log re), 



(5ire 

62 log re log re 


-1 


and g^nfi^Rn^ei) = g'^S.^ ^(logre — (52loglogre). Using Proposition 2.6 (i) in 


log 




(5ire 


log re — (52 log log re 


log 


(5ire 


log re — (52 log re log re 


1 -1 


da ^ — g'^S^ ^ < 0 , 


dx-l 


re 


00 . 


7'^(5 


-1 

1 


At the last inequality, we applied the constraint in (5i. Therefore, we have dlog Rn'^ — g^nf (^Rn'^ ei) 
— 00 , re —>■ 00 , as requested. □ 
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Proof, {proof of Proposition \ 6 .S^ Since G RVii^., it is easy to show (I6.10p . and therefore, we 
prove only that dlog-Rn^ — g'^nf{Rn^ei) —>■ —oo as n —>■ oo. We see that 

d log R^f^ < d log (log n) ~ dv~^ log log n , n —>■ oo , 

and that g^nf {Rn^eij = loglogn. By virtue of the constraints in hi and 62 , we have 

dv~^ — < 0; thus, the claim is proved. □ 

Remark 6.4. The proof of Lemma 16.31 merely estimated the probability in (16.lip from below. 
Therefore, it seems to be possible that in the propositions above, (16.411 and (16.9p may hold for the 
sequence Rn 00 growing more quickly than Rn^ but more slowly than Rlf^\ i.e., R^f^ < Rn ^ 
R^'^] it is unknown, however, to what extent we can make Rn closer to 


7. Proof of Main Results 


This section presents the proof of the main results of this paper. The proof is, however, rather 
long, and therefore, it is divided into several parts. All the supplemental ingredients necessary are 
collected in the Appendix, most of which are cited from [211- 

Let Ann(iL, L) be an annulus of inner radius K and outer radius L. For xi,..., Xfc G M'^, define 
Max(xi,..., Xfc) as the function selecting an element with the largest distance from the origin. 
That is, Max(xi,... ,Xk) = Xi if |lxj|| = maxi<j<fc ||xj||. If multiple Xj’s achieve the maximum, we 
choose an element with the smallest subscript. 

In the following, T, etc. always represent a finite collection of d-dimensional real vectors. 

We use the following shorthand notations. That is, for x = (xi,...,Xm) G x G and 

y = (yi,...,y^_i)G(M''r-i, 

/(x) := f{xi) ■ ■ ■ f{Xm) , 
f{x + y) := f{x + yi)--- f{x + ym-i) , 
h{ 0 ,y) := h{ 0 ,yi,... ,ym-i) etc. 


Regarding the indicator ht : —>■ {0,1} given in (j2.ip . the following notations are used to 

save space. 


(7.1) 


:= hf (x) — /is(x), 0 < s <t, X G (M'^)^, 

:= hf (x) - /i^(x), 0 < s < t, X G (M'^)^, 

:= ht^s{x.) l{m(x) > R„} , 0 < s < t, x G (M'^)*', 


ht,s(x) := 

:= 

(7.2) Vm(x) := 
and for £ G {0,..., k}, 

(7.3) h[^l{yL) := ht{xi,... ,Xk)hs{xi,... ,Xi,Xk+i,... ,X 2 k-£), Ls>0, x G 
In particular, we set 

, . ^ ] hs{xk+i,.. ■ ,X2k) if^ = 0, 

• • • ? • • • ? ^2k—i) •— \ 1 / \ -c n 1 

\^1 5 • • • 7 ^k ) 11 -c — k . 

In Section 17711 we use, for 1 < iL < L < 00, n G N+ and t > 0, 

hn,t,K,Li^) ■= ht{x) l{m(x) > Rn , Max(x) G Ann{KRn, LRn)} , 
hn,t,K,Li^) •= W l{”i(x) > Rn , Max(x) G Ann{KRn, LRn)} . 
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The same notations are retained for Section O to represent, for 0 < K < L < oo, n G N+ and 
t > 0, 

hn,t,K,L{x.) := ht{x) l{m(x) > Rn, a{Rn)~^ {Max{x) - Rn) G [K,L)] , 

•= > Rn , (Max(x) - Rn) G [K,L)} . 

Finally, C* denotes a generic positive constant, which may change between lines and does not 
depend on n. 

In the following, we divide the argument into two subsections. Section ITT] treats the case in which 
the underlying density has a regularly varying tail; our goal is to prove Theorem 14.31 Subsequently 
Section[72]provides the proof of Theorem l5.3[ where the density is assumed to have an exponentially 
decaying tail. Before the specific subsections, however, we show some preliminary results, which 
are commonly used in both subsections for the tightness proof. 

Lemma 7.1. Let ht : —)■ {0,1} be an indieator given in (|2.1I) . Fix T > 0. Then, we have for 

£ € {1,..., k}, 

(7.4) f dy [ dz 2 [ y, zi)/i+^(0, y, Z2) < - s)(s - r), 

[ dy [ dz 2 [ dzi/i7j0,y,zi)/i“^(0 ,y,Z 2 ) < C*(t-s)(s-r) 

./(R‘')^-i 4(R'*)'=-< 

for aU0<r<s<t<T. 


Proof. We only prove the first inequality. 11 £ = 1 or £ = k, the claim is trivial, and therefore, we 
can take 2 < £ < k — 1. It follows from (j2.6p that the integral in (j7.4p is not altered if the integral 
domain is restricted to [B{0, kT))^ ^ x [B{0, kT))^ ^ x (B(0, kT))^ ^. With A being the Lebesgue 
measure on we see that for every y G 

(7.5) [ (0,y,z)(iz = A{z G (B(0,A;r))^“^ :/i)^(0,y,z) = 1, /i+(0, y, z) = O} 

J[B{ 0 ,kT)) 

< A{z G (B(0, kT))^ ^ s < \\zi — ZjW < t lor some i ^ j) 

+ A{z G (B(0, kT))^ ^ : s < Hzj — yj|| < t for some i,j) 

+ A{z G (B(0, kT))^ ^ : s < HziH < t for some i] 

+ A{z G [B{0,kT))^ ^ : s < \\yi — yj\\ < t for some i ^ j} 

+ A{z G (B(0, kT))^ ^ : s < ||yj|| < t for some i]. 

Observe that for i ^ j, 

A{z G (B(0,fcr))^"^ : s < \\zi - ZjW < t} < - s ^), 

where iVd is the volume of the d-dimensional unit ball. Since the second and the third terms on the 
rightmost term in (17.51) have the same upper bound, we ultimately obtain 


'(B(0,fcT))' 


f^M(0:y:Z)dz 

. i-l i-l 

<C*{t'^-s‘^+ ^ l{s < ||yi-yjll < t}+ ^l{s < llyill < t} ). 

*=i 
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c 


I 


{B{0,kT)) 


Therefore, the integral in (17.41) is bounded above by 

. e-i £-1 

l{s < ||2/i< t}+ ^l{s < llyill < t} 

, i-i i-\ . 

X Y < lli/i-yjll < s}+ 

^ i=l '' 

An elementary calculation shows that for all G {1,... — 1} with i > j and i' > f, 

L st-i ij'S <\\yi- VjW < t} l{r < Wvi^ - yfW < s}dy 

J[B{0,kT)) 

< C*{t'^ - s^){s'^ - /) < C*{t - s){s - r) 

In particular, if i = i' and j = j', the integral is identically zero. Applying the same manipulation 
to the integral of other cross-terms, we can conclude the claim of the lemma. □ 


7.1. Regularly Varying Tail Case. Under the setup of Theorem l4.,'l[ we first define the subgraph 
counting process with restricted domain. For 1 < it' < L < oo, n G N+, and t > 0, let 


and 


Gn,KAt) = Y ^*0^) l{m(T) > Rn , Max(T) G kTin{KRn,LRn)} 
yc.v„ 

■= Y^ hn,t,K,L{y) , 

ycVr, 


^{m{y) > Rn , Max(T) G Ann{KRn,LR^)} 

ycVn 

■■= E . 

ycVn 

where (Rn) satisfies (14.2p . For the special case K = I and L = oo, we simply denote G„(t) = 
Gn,i,oo{t) and G^(t) = G^^ oo(^)- subgraph counting processes, centered and scaled, for which 
we prove the FCLT, are given by 

Xn{t)=r-^G(^Gn{t)-E{Gn{t)}), 

(7.6) Xtit) = rY^Gtit) - lE{G±(t)}) , 

where (r^) is determined by (14.3p according to which regime is considered. The first proposition 
below computes the covariances of {Gn,K,L{t)) ■ 


Proposition 7.2. Assume the conditions of Theorem f.3. Let 1 < K < L < oo. 

(i) If nfiyRnCi) —)• 0 as n ^ oo, then 

r-i Cov{Gn,K,L{t), Gn,K,L{s)) ^ ” L‘^-‘^^)Lk{t, s), n^oo. 

(a) If nf{Rnei) —^ G (0, oo) as n ^ oo, then 

T-^Cov{Gn,KAt)^Gn,KAs)) ^ ^_ L^-^^^^-^'>)f’^-^L,{t,s) , 


n ^ OO. 
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{iii) If nf{Rnei) —>■ oo as n —>■ oo, then 

T~^ Cov{Gn,K,L{t), Gn,K,L{s)) ^ n ^ OO . 

Proof. We start by writing 

^{Gn,K,L{t) Gn,K,L{s)} 
k 

= EK E E l(3^i)/ in,., r,l( 3^2) 1{ 13^1 n 3^21 = ^ }} 

£=o yicPny2cVu 
k 

■.= J2nh}. 

£=0 

For i = 0, applying Palm theory (see the Appendix) twice, 

jf^k 

E{3o} = J^^^{hn,t,K,LiXl, ■ ■ ■ ,Xk) hn,s,K,L{Xk+l, ■ ■ ■ ) A'2fe)} 

= ^{Gn,K,L{'t)}^{Gn,K,L{s)} ■ 

Therefore, the multiple applications of Palm theory yield 

k 

Cov[Gn,K,L{t) ,Gn,K,L{s)) = 


£=1 


k ^2k-£ 


E 


■E{hn,,,K,L{yi)f‘n,>.K,L{y2) R iJl nJ'jl = <}}. 


fcr <!((«:-<)!)" 

Define for ^ € {1,... , k}, 

G^J:\K, L) := {x G : Max(xi,..., Xfc) G Ann(A:i?„, LRn) , 

Max(xi,... ,X£,a;fc+i,... ,X 2 k-£) G Ann{KRn, LRn)} . 

By the change of variables x —>■ (xjX + y) with x G x G y G together 

with invariance (j2.2p . while recalling notation (17.31) . 

E<! 


X 


'\hn,t,K,L{yi) hn,s,K,L{y2) l{ |3^1 H T 2 I = ^ } } 

/(x) l{m(x) > Rn} /i|2(x) l{x G C'W(A:,L)}d 

= / / f{x)f{x + y)l{m{x,x + y)>Rn}h[^l{0,y) 

X l{(x,x + y) G G}^\K,L)}dydx. 

The polar coordinate transform x ix,0) and an additional change of variable p —>■ r/Rn yield 

(7.7) E{Vi,i^,L(3^i) VAifT(3^2)i{ 13^1 nTsI =^}} 

= RifiRnei)^’^-^ [ J{9)d9 r dp [ 

JSd-i 3i 3(Rd)2fe-i-i j[Rnei} 

n l{\\p9 + y^/Rn\\ > l}/iig(0,y) 


2=1 


X 
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X l{{Rnpe,R^pe + y)eC^\K,L)], 


where Sd-i denotes the {d — l)-dimensional unit sphere in and J{9) is the usual Jacobian 


J(6') =sin^ ^(6»i)sin^ ^(6*2) • • • sin(6'fc_2). 


Note that by the regular variation of / (with exponent —a), for every p > 1, 9 £ Sd-i, and y^’s, 


(7.8) 


fjRnpei) 

fiRnei) 


TT f{\\RnP9 + UiWei) a(2k-£-l) 


n —)■ 00 


and, furthermore, 

(7.9) l{iRnp9,Rnp9 + y) £Ci^\K,L)} ^1{K <p<L}, n ^ 00 . 


Substituting (j7.8p and (17.91) back into (j7.7n . while supposing temporarily that the dominated con¬ 
vergence theorem is applicable, we may conclude that 


(7.10) Cov{Gn,K,L{t) ,Gn,K,L{s)) 

k 

s), n ^ 00 . 

£=1 


Observe that the limit value of nf{Rnei) completely determines which term on the right hand side 
of (j7.10p is dominant. If nf{Rnei) —)• 0, then the fcth term,i.e., i = k,m the sum grows fastest, while 
the first term, i.e., i = 1, grows fastest when nf{Rnei) —)■ 00 . Moreover, if nf{Rnei) —^ G (0, 00 ), 
then all the terms in the sum grow at the same rate. This concludes the claim of the proposition. 

It now remains to establish an integrable upper bound for the application of the dominated 
convergence theorem. First, condition ()2.3p provides 

/lifi (0,y) < l{||yi|l <kit + s), i = l,-..,2k-i-l} . 

Next, appealing to Potter’s bound ,e.g.. Proposition 2.6 {ii) in [23], for every ^ € {0,a — d) and 
sufficiently large n. 

Up > 1} < (1 + {) p-“+« Up > 1) 

f{Rnei) 

and 

“ii"' + 1{HP0 + ydRM > 1} < (1 . 


Since p^ ^ < 00, we are allowed to apply the dominated convergence theorem. □ 

The next proposition proves the weak convergence of Theorem 14.31 in a finite-dimensional sense. 

Proposition 7.3. Assume the conditions of Theorem \4-S\ Then, weak convergences (i) — (Hi) in 
the theorem hold in a finite-dimensional sense. Furthermore, let be the processes defined in 
(HSj). Then, the following results also hold in a finite-dimensional sense. 

{i) If nf{Rnei) —>-0 as n ^ 00 , then 

(7.11) 


(X+,X“)^(V+,V-). 
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(m) If nf{Rnei) f, € ( 0 , oo) as n ^ oo, then 


(7.12) (X+,X-) ^ ■ 

\£=1 l=l / 

{Hi) If nf{Rnei) —>■ oo as n ^ oo, then 

(7.13) (X+,X-)^(V+,Vr). 

The limiting Gaussian processes (V^, V^), I = 1,... ,k are all formulated in Section \472 . 

Proof. The proofs of ()7.1ip . (I7.12p . and (I7.13P are a bit more technical, but are very similar to the 
corresponding results in Theorem l4.31 therefore, we check only finite-dimensional weak convergences 
in Theorem 14.,31 The argument here is closely related to that in Theorem 3.9 of |21] . for which we 
rely on the so-called Cramer-Wold device. For 0 < < ■ • • < < oOj Rij ■ ■ ■) € M and m > 1, 

define Sn '■= cijGn{tj). For K > 1, Sn can be further decomposed into two parts: 

m m 

Sn — ^ ^jGn,l,K {tj) T (^jGn,K,oo{ij) 
i=i i=i 

• = Tji^) 

n ' 

We define a constant 'Jk as follows in accordance with the limit of n/(ii„ei). 

ifra/(72„ei) ^ 0 , 

IK := ET=1 Er=i Eii(l - t,) if nf{Rne,) ^ ^ G ( 0 , oo), 

Z)jLi ®i«i(l - ifn/(i 2 nei) ^ oo . 

Moreover, 7 := limi^_^oo 7 ^^. It is then elementary to check that, regardless of the regime we 
consider, 

^Var{TW}yA', r"^Var{17 ,(^)}7 - 7 i^ as n 00 . 

For the completion of the proof, we ultimately need to show that 

r-i/2(5„_E{54) ^iV(0,7). 

By the standard approximation argument given on p. 64 of m, it suffices to show that 

(7.14) 


^ N( 0 , 7 i^) for every K> 1 ; 


equivalently, 

( 7 . 15 ) 


rW_E{Tr^} 


(K)^ 


A^( 0 , 1 ) for every K > 1 . 


Let (^Qi • ^ G 


^Var{Ti^)} 

be a collection of unit cubes covering W^. Define 
:= {£ G N : n Ann(i?„, KRn) / 0} , 

where we have that | 14 | < C*Rf,. 

Then, Tn^^ can be partitioned as follows. 

m 

E E% E (T) l{m{y) > Rn, Max(T) G Ann{Rn,KRn) n Q^} 

teUn i = l yCPn 
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:= Y1 • 

teVn 

For i,j G Vn, we put an edge between i and j (write z ~ j) if z 7 ^ j and the distance between Qi 
and Qj are less than 2ktm- Then, gives a dependency graph with respect to ^ £ 14); 

that is, for any two disjoint subsets Ii, I 2 of Vn with no edges connecting Ii and I 2 , {rji^n, ^ £ h) 
is independent of {rji^n, ^ £ h)- Notice that the maximum degree of iVn,^) is at most finite. 

According to Stein’s method for normal approximation (see Theorem 2.4 in [21]), (|7.15|i imme¬ 
diately follows if we can show that for p = 3,4, 


(7.16) 


Rt max 




0 as n —>■ 00 . 


Since the proof for showing this varies depending on the limit of nf^RnCi), we divide the argument 
into three different cases. Suppose first that nf{Rnei) —)■ 0 as n —>■ 00 . Let Zi^n denote the number 
of points in Vn lying in 

Tube(Q£; ktm) '■= {x G : inf ||x — y|| < ktm} ■ 

y^Qt 

Then, Zg^^n has a Poisson distribution with mean n frpnhe(Qe-kt ) Using Potter’s bound, 

we see that Zg^n is stochastically dominated by another Poisson random variable Zn with mean 
C*nf{Rnei). Observing that 


l%„l < c* (^'•") . 

we have, for g = 1,2, 3,4, 

^r]g,nV < U*E < C*E 


<C*{nfiRnei))\ 


where in the last step we used the assumption nf{Rnei) —>■ 0 . 
It now follows that for p = 3,4, 


Therefore, 


max£'|ry£_„ - E{r]g^n}f < C*{nf{Rnei)) 

Vn 






(varjri^)}) 


p/2 


{n^Rif{RneiYlK) 


p/2 


c 


— {n^Rif{R^egf) 
7k 


kV-pV 


0 . 


n 


00 , 


where the last convergence follows from (|4.2p . 

In the case of nf{Rnei) ^ ^ G ( 0 , 00), the argument for proving (I7.16h is very similar to, or even 
easier than, the previous case, so we omit it. 

Finally, suppose that n/(i?„ei) ^ 00 as n —>■ 00. We begin by establishing an appropriate upper 
bound for the fourth moment expectation 

4 


(7.17) 


- E{rjg,n}\ = E ( / ) (-l)^E{E}(E{%,n}) 

j=0 


4-i 
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Letting 

m 

9^,n{y) •= l{m(3^) > Rn-, Max(3^) G Ann(ii:„, n Qij , 

i=i 

we see that for every j € {0,..., 4}, 

Fn{j) :=E{ry^;j(E{%, 4 )"-^' 
can be denoted as the expectation of a quadruple sum 


(7.18) E E E E E 9e,n{yi) 9e,n{y2) gi,n{y3) ge,n{y4) z , 

I ^1 y2yA J 

where each of Vn \ ■ ■ ■ ,'Pn'^ is either equal to or an independent copy of one of the others. By 
definition, each 3^j is a finite collection of d-dimensional vectors. If, in particular, |3^iU3^2U3^3U3^4| = 
Ak, i.e., any two of 3^*, i = 1,... ,4 have no common elements, then the Palm theory given in the 
Appendix reveals that (j7.18j) is equal to (E{? 7 £^„}) . Then, in this case, their overall contribution 
to ()7.17p is identically zero, because 

j=o 


Next, suppose that iTi U T 2 U Ts U 3^41 =4/0 — 1, i.e., there is a pair (Tj, Tj), i ^ j having exactly 
one element in common and no other common elements between Tj’s are present. In this case, 
(j7.18l) can be written as 


(7.19) 


^ 2 /c 1 ^ 

--^ E ge,n{yi) 9e,n{y2) l{ |3^i n T 2 I 

((/c- 1 )!) 



In particular, ()7.19l) appears once in Fn{2), 
total contribution to (j7.17l) sums up to 


times in i3i(3), and 


times in i3i(4). Thus, the 







X (17dm = 0 . 


We may assume, therefore, that |TiUT2U3^3UT4| < 4/c —2. Let us start with |TiUT2U3^3UT4| = 
4/0 — 2, where we shall examine in particular the case in which Vn^ = Vn^ = = 'P^\ 

|3^i n T 2 I = 2 and no other common elements between Tj’s exist. The argument for the other cases 
will be omitted because they can be handled in the same manner. Then, by Palm theory, (I7.18|) is 
equal to 


(7.20) 


n 


2 k-2 


2 ((A:- 2 )!) 


■E{dr,n(3^i)5£,n(3^2)l{|3^inT2| = 2 }} (^E{ff,,„(T)} 


Because of Potter’s bound, together with the fact that Qi intersects with Ann(i?„, iLi?„), 
E{ff^,n(h3)ff£,n(3^2)l{|3^1 nT 2 | = 2 }} 

< C*(p{Xi G Tube(Q,;A:t^)})'''”' < C*f {Rnei)^^-^ 
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Similarly, we can obtain 


^geAy)} <C*f{RneA 


and therefore, the absolnte value of ()7.2np . equivalently that of ()7.18p . is bounded above by 

C*{nf{Rnei)y^-\ 

A similar argument proves that if |3^i U 3^2 U 3^3 U 3^4| = 4k — q for some g > 3, the absolute value 
of (|7.18l) is bounded above by C* [nf{Rnei))'^^ Putting these facts altogether, while recalling 
nf{Rnei) ^ oo as n —>■ oo, we may conclude that 


E| 77 £,„ - < C*{nf{Rnei))^'^ \ 

Now, it is easy to check (I7.16p . 

In terms of the third moment expectation — E{ry£^„}|^, we apply Holder’s inequality to 

obtain 

EA,n - E{rje,n}f < < C*{nf{Rnei)f’'-^^\ 

Again, it is easy to prove ^7.161) . 

Now, we have obtained a CLT in (I7.14p as required, regardless of the limit of n/(i?„ei). □ 


An important claim is that once the tightness of each and X“ is established in the space 
P[0, oo) which is equipped with the Skorohod Ji-topology, the proof of Theorem 14.31 is complete. 
To see this, suppose that X+ and X“ were both tight in P[0, oo). Then, a joint process (X+, X“) 
is tight as well in P[0, oo) x P[0, oo), which is endowed with the product topology. Because of the 
already established finite-dimensional weak convergence of (X+, X“), every subsequential limit of 
(x+,x- ) coincides with the limiting process in Proposition 17.31 This in turn implies the weak 
convergence of (X+,X“) in P[0,oo) x 2?[0, oo). Using the basic fact that the map {x,y) ^ x — y 
from P[0, oo) x P[0, oo) to P[0, oo) is continuous at (x, y) € C[0, oo) x C[0, oo), while recalling that 
the limits in Proposition 17.31 all have continuous sample paths, the continuous mapping theorem 
gives weak convergence of X„ = X+ — X“ in T>[0, oo). 

Proposition 7.4. The sequences (Xj)^) and (X“) are both tight mP[0,oo), irrespective of the limit 
of nf{Rnei). 


Proof. We prove the tightness of (X+) only, in the space P[0, L] for any fixed L > 0. For notational 
ease, however, we omit the superscript “-I-” from all the functions and objects during the proof. 
By Theorem 13.5 of [8], it is sufficient to show that there exists B > 0 such that 

E{(A4f) - Xr,{s))\Xn{s) - Xnir))^] < B{t - rf 

for all 0 < r < s < f < L and n > 1. 

For typographical convenience, we use shorthand notations (17.11) . (17.21) . and further. 


Then, 

E{(X„(f) 




VmW- 

ycVr. 


Xn{s)f{Xn{s)-Xn{r)f] = 




E{en.s,r-})'} 



{-lY+<iF^{p,q), 
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where 


Note that for every p,q € {0,1, 2}, Fn{p, q) can be represented by 


(7.21) 


E 


E E E E hn,t,siyi) hn,t,s{y2) hn,s,riy 3 ) hn,s,r{y 4 :) > , 


y2cpi^'’ yacvL^'^ y^cvi*'* 


where each of Vn ^, ■ ■ ■, is either equal to or an independent copy of one of the others. 

According to the Palm theory given in the Appendix, if |3^i U 3^2 U 3^3 U 3 ^ 4 1 = Ak, i.e., any two 
of X have no common elements, then (17.211) reduces to Then, an overall 

contribution in this case identically vanishes, since 

2 2 


EE 

p =0 q =0 


(_l)P+a(E{^„, = 0. 


In the following, we examine the case in which at least one common element exists between Ti’s. 
First, for i = 1,... ,k, we count the number of times 

(7.22) E (E E'>»a s{yi) K,tAy2) 1 { ITi n T2I = ^}} (E{en,.,.})' 

yiCVn y2CVn 

appears in each Fn{p,q). Indeed, (I7.22h appears only once in F„(2,0), F„(2,1), and F„(2,2). 
Therefore, the total contribution amounts to 




(_l)2+0 ^ 




(_1)2+1 + 




X (n:m = o. 


Similarly, for every i = 1,... ,k, no contribution is made by 

E E r-(Tl) K,sAy2) 1{ |Tl n T2I = ^}} (E{enA.})'. 

yiCVn y2CVn 


Subsequently, for ^ = 1,..., /c, we explore the presence of 

(7.23) E { E E s(Tl) V.,r-(T 2 ) l{|Tl n T 2 I = • 

yiCPn y2CVn 

One can immediately check that (17.231) appears once in F„(l, 1), twice in F„(2,1), twice in Th(l, 2), 
and four times in F„(2,2). However, their total contribution disappears again, because 


(- 1 ) 


1+1 


+ 


(_l)2+i . 2 


+ 


(- 1 ) 


1+2 


2 + 


(_l)2+2.4 


X (Tt:^ = 0 , 


Next, let G {0,..., k}, i = 1, 2,3, G {2,... , 2k} such that at least two of ii's are non-zero, so 
that we should examine the appearance of 


(7.24) E (E E E ^n,t,s(Tl) hn,t,s{y2) hn,s,r{y3) 

yiCVn y2CVn ysCVn 

X i{|Ti n T2I = h, |Ti n TsI = ^2, IT2 n T3I = A, |Ti u T2 u TsI = 3 * - £}} . 
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This actually appears once in Fn{2, 1) and twice in Fn{2, 2); therefore, their overall contribution is 


1 (- 1 ) 


2+1 


+ 


(_l)2+2 . 2 


X ^nm = 0, 


For the same reason, we can ignore the presence of 


E E hn,t,s{yi) hn,s,r{y 2 ) hn,s,r{y3) 

yiCVn y2CVn ysC^Vn 

X i{|Ti n T2I = h, 13^1 n TsI = h, 13^2 n T3I = h, | 3 ^i u T2 u T3I = ^k- £}} E{en,t,.}. 

where G {0,..., k}, i = 1,2, 3, ^ G {2,..., 2k} such that at least two of £j’s are non-zero. 

Putting these calculations altogether, we find that the tightness follows, once we can show that 
there exists B > 0 such that 


(7.25) r-^Ej E E E E hn,t,s{yi) hn,t,s{y 2 ) ^n,s,r(3^3) ^n,s,r(3^4) 

yiCVn y2CPn ysCPn y4.CV„ 

X 11 each 3’j has at least one common elements with 

at least one of the other three} | < B(t — r)^ 

for allO<r<s<t<L and n > 1. We need to check only the following possibilities. 

[I] ^ := |Ti n T2I G {1,..., k], i' := IT3 n T4I € {1,... , k}, and (Ti U T2) n (T3 U T4) = 0- 
[D] £:= |T2nT3| G {!,..., A:}, := |TinT4| G {!,...,A:}, and (T2 U T3) n (Ti U T4) = 0- 

[in]. Each T’j has at least one common element with at least one of the other three, but neither [I] 
or [n] is true. 

For example, if |3^i n 3’2| = 2, |3’i n 3’3| = 3, [3^2 H 3’4| = 1, and there are no other common 
elements between 3^i’s, then it falls into category [m], where, unlike [I] or [II], the expectation in 
(|7.25p can no longer be separated by the Palm theory. 

Denoting by A the left-hand side of ()7.25p , let us start with case [I]. As a result of Palm theory, 

2k—I 

^ Vm(3^i) K,tAy2) i{ I3^i n T 2 I = 

l\[\k — i)\j '' '' 

9h—f' 

X T-^ _ , E[hn,s,rm K^sAy^) lIl^S H T 4 I = /}} 

I i j!j 

:= Ai X A 2 . 


Proceeding as in the calculation of Proposition 17.21 we obtain 

(7.26) [ dy [ dz 2 f dzi ht,^(0, y, zi) ht,^(0, y, Z 2 ), 

J{K.d-)k-e J(]Rd)fe-£ 

{1.27) A 2 <C*T~^n^'^~^'Rif {Rneif'^~^' j dy[ dz 2 [ dzi/is_^(0,y, zi)/is_^(0, y, Z 2 ). 

J(Udy'-i 

Notice that ht is increasing in t in the sense of (12.51) (recall that the superscript “-I-” is suppressed 
during the proof). It also follows from (12.611 that the triple integral in (17.261) is unchanged if the 
integral domain is restricted to {B{0, kL))^ ^ x (^B{0, kL))^ ^ x {B{0,kL))^ ^. Therefore, with A 
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being the Lebesgue measure on 

/ dy dz2 / dziht,s{0,y,zi) ht^s{0,y,Z2) 

< X{{B{0,kL))^~^} f [ ht,s{0,y,z)dydz 

J{Rd.y-'^ JiRd^k-e 

= X{ (5(0, kL)) [ h{0, y)dy 

J(]Rd)fc-l 

< C*{t-r). 

Applying the same manipulation to the triple integral in (I7.27p . we obtain 

A < - rf. 

It remains to check that sup„ R"^f{Rnei)‘^^~^~^' < oo, which is, however, easy to prove, 

irrespective of the definition of Tn- Now case [I] is done. 

Next, we turn to case [It]. As a consequence of the same operation as in [I], we obtain the same 
upper bound for A up to multiplicative constants. 

Finally, we proceed to case [IH]. Let i := Ak— |3^i U 3 ^ 2 U 3 ^ 3 U 3 ^ 4 l; then, it must be that 3 < i < 3k. 
It follows from Palm theory that 

A = c*T-^n^^-^E{K,t,s(yi) K,tAy2) hn,sAy3) KsAy^)} 

with (3^1,..., 3 ^ 4 ) satisfying requirements in case [IH]. In particular, ( 3^1 U 3 ^ 2 ) H ( 3^3 U 3 ^ 4 ) must be 
non-empty; hence, we may assume without loss of generality that 3^i n 3^3 / 0. Set i' := |3^i n 3 ^ 3 | G 
{!,..., k}. By (12.51) and (12.61) . we have 

A<C*TAA^~^Rif{Rneif^~^ [ dy [ dz 2 [ dziy, zi)y, Z 2 ), 

J(Rd.y'-l J(Rd)fc-£' 

Because of Lemma 17.11 

A < C*TAA^-^Rif{Rneif^-\t - rf. 

Once again, verifying 

supr-2n^^-^5^/(5„ei)^^-^ < 00 

n 

is elementary, and hence, we have completed the proof of (I7.25P as required. □ 

7.2. Exponentially Decaying Tail Case. We start by defining a subgraph counting process with 
restricted domain. For 0 < TL < L < 00 , we define 

Gn,KAt) = Y. ht{y)Arn{y) > 5, , a(5„)-i (Max(T) - 5,) G [K,L)] 
ycVr, 

■= hn,t,K,L{y) , 

yCVn 

Arn{y) > Rn, o(5„)-i(Max(T) - R^) G [K,L)} 

ycVn 

■■= Y hYxAy ^, 

ycVn 


and 
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where (Rn) satisfies (j5.7l) . For the special case K = 0 and L = oo, we denote Gn{t) = Gn,o,cxi{t) 
and G^{t) = G^q ^(t). The centered and scaled versions of the subgraph counting process are 

(7.28) Xnit) = T-^/^(Gnit) - lE{G„(t)}) , 

(7.29) X±(t) = r-V2(^G±(t) - ]E{G±(t)}) , 

where (r^) is given in (15.8p . As seen in the regularly varying tail case, we first need to know 
the growing rate of the covariances of Gn,K,L{t)- Before presenting the results, we introduce for 

s) := Di r [ E£r^-^(ei,.d 

Jo 

X l{y G E^^]^{p,ei)} h[^l{0,y) dydp, t,s>0, 
where Di is given in (15.91) . h[^l {0,y) is defined in (15.111) . and for p > 0 and 6 G Sd-i, 

^) = {y € :p + c-\9,y,)>0, i = l,...,2k-i-l, 

A'< maxjp, p + 

K < max(p, p + c~^ max (9, Vi)] < L > . 

Note that M^ o,oo(A'S) completely matches (15.101) . 

Proposition 7.5. Assume the conditions of Theorem \5.,‘A Let 0 < K < L < oo. 

(i) If nf{Rnei) —)■ 0 as n ^ oo, then 

T~^ Cov[Gn,K,Lit), Gn,K,L{s)) Mk,K,L{t, s) , n ^ OO . 

(a) If nf{Rnei) —^ G (0, oo) as n ^ oo, then 

k 

Cov{Gn,K,L{t), Gn,K,L{s)) ^ ^ s) , U ^ OO . 

£=1 


{Hi) If nf{Rnei) oo as n ^ oo, then 

T~^ Cov{Gn,K,L{t), Gn,K,L{s)) s) , U ^ OO . 

Proof. As argued in Proposition 17.21 with the multiple applications of Palm theory, one can write 


2k-e 


Cov{Gn,K,L{t) ,Gn,K,L{s)) =y^- , 
Define for ^ G {1,... , k}, 


^{hn,t,K,L{yi) K,s,K,L{y2) l{ |Tl H Ts | = ^ } }. 


F^^\K,L) := {xG (M'')2^-^a(A„)-i(Max(xi,...,Xfc)-A„) G [K,L), 

a(An)“^(Max(xi, .. .,xi,Xk+i, ■. ■,X2k-e) - Rn) G [K,L)} . 

By the change of variables x ^ {x,x + y) with x G x G y G (]R‘^)^^“^“^, together 

with invariance (j2.2p . 

^{hn,t,K,L{yi) hn,s,KAy2) 1 { |Tl H ^21 = ^ } } 
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/ /(x) l{m(x) > Rn} hfl{yi) l{x e FW(K,L)}dx 

/ / /(a;)/(x + y)l{m(a:,x + y) > iZ„}/ij5(0,y) 

X l{(a;,a; + y) G F^\k, L)]dydx . 


Let Jk denote the last integral. Further calculation by the polar coordinate transform x —>■ (r, 9) 
with J{9) = \dx/d9\ and the change of variable p = a{Rn)~^{r — Rn) yields 

r 1*00 p 

(7.30) Jk =a{Rn)Rt^f{Rnei)^^-^ / J(0)d0 dp dy 

Jo ./(Rd)2''-^-i 

/((^n + a(7?n)p)ei) 

2fc-£-l 

X f{Rnei)~^f(^\\{Rn + a{Rn)p)9 + yi\\ei'^l^\\{Rn + a{Rn)p)6 + yi\\>Rn^ 

i=l 

xl{{{Rn + a{Rn)p)e, {Rn + a{Rn)p)e + y) GFW(iL,L)}hg(0,y), 
where S'^~^ is the {d — l)-dimensional unit sphere in 

The following expansion is applied frequently in the following. For each i = 1,... ,2k — £ — 1, 


{Rn + o,{Rn)p)0 + yi 


Rn + a{Rn)p + {9, yi) + 7„(p, 9, y^), 


so that 7n(/0, 9, yi) —)• 0 uniformly in p > 0, 0 G 5'^“^, and ||?/j|| < k{t + s). 

For the application of the dominated convergence theorem, we need to compute the limit of 
the expression under the integral sign, while establishing an integrable upper bound. We first 
calculate the limit of the indicator functions. For every p > 0, 0 G , and ||yil| < k{t + s), 
i = l,...,2k-£-1, 


2k-l-l 

'^^\{Rn 0,{Rn)p)9 + yi\\ > Rri^ 

i=l 

xl[{{Rn + a{Rn)p)9, {Rn + a{Rn)p)9 + y) GTW(i^,L)} 
^ ^ > n^oo. 


Next, it is clear that for every p > 0, (l + a{Rn)p/Rn)^ ^ tends to 1 as n —)• oo (see (j5.3p i and 
is bounded above by 2(max{l,p})^^ 

As for the ratio of the densities in the second line of (j7.30l) , we use the basic fact that 1 /o is flat 
for a, that is, as n —>■ oo, 


(7.31) 


(k{Rn) 

a{Rn + a{Rn)v) 


1) 


uniformly on bounded ?;-sets; 


see pl42 in m for details. Noting that L is also flat for o, we have for every p > 0, 


/((/!„+ o(fl.,)p)ei) L { Rr . + a { Rr .) p ) f , ,o^^, ,.bi 1 

-7(17i)-“- WQ - ^■S<[-i>(R..+a{R^)p)+iiRn)) 
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_ L[Rn + a{Rn)p) f a{Rn) , I 

L{Rn) Jo a{Rn + a{Rn)r) i 

e~^, as n —>■ oo . 


To provide an upper bound for the ratio of the densities, let {qm{n), m > 0,n > l) be a sequence 
dehned by 

qm{n) = a{Rn)~^ {i^{Rn) + m) - RrJj , 

equivalently, 

ijj{Rn + a{Rn)qm{n)) = ipiRn) + m . 

Then, for e G (O, (d + 7(2A: — ^))“^), there exists an integer A^e > 1 such that 

qm{n) < for all n > iVg, m > 0 . 

For the proof of this assertion, the reader may refer to Lemma 5.2 in [3]; see also Lemma 4.7 of 
|20j . Because of the fact that ijj is non-decreasing, we have, for sufficiently large n, 

exp + a{Rn)p) + V'(-Rn)| l{/0 > 0} 

oo 

= ^ l{qmin) < p< qm+iin)} exp|-'0(i?„ -h a{Rn)p) + V'(-Rn)| 
m=0 
oo 

< ^ 1{0 < p < e-^. 

m=0 


Using the bound in (j5.5p . 


L{Rn) ^L{Rn +a{Rn)p)l{p > 0} < C 



7 

< 2C'(max{p, 1})"*^. 


Combining these bounds. 


f(^{Rn + a{Rn)p)ei^ 
fiRnCl) 


l{p > 0} < 2C(max{p, 1})^ ^ l{0 < p < e"™ 

m=0 


Finally, we turn to 


2k-i- 


/(iKflr. + a(Rn)p)e + wild) 2M-1 L(fl„ + a(R„){p + („{p, 8, w))) 


n 

2=1 


/(A.ei) 


n 

2=1 


X exp < — 




L{Rn) 

P+iApfi,yi) q(^^) 


O^i^Rn T ^i,Rri)^^ 


dr } , 


where 


Since c = hmn_>.oo a{Rn) > 0, 


^n{p,0,y) = 


(<9,p) + 7n(p,6>,p) 

(ki^Rn) 


A := sup \Cn{p,d,y)\ < oo . 

n>l, p>0, 

0eS‘^-\ |L||<fc(t+s) 
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Therefore, because of the uniform convergence in (I7.31|) . for every p > 0, 6 £ S'^ and ||yj|| < 
k{t + s), 


2k-i- 


n 


^-1 f(^\\{Rn + a{Rn)p)0 + ViWei^ 


2k-l-l 


f{Rnei) 


exp{-{2k - i-l)p-c ^ ^ {0,yi)}. 


i=l ^ ' i=l 

Subsequently, on the set 

^\\{Rn + a{Rn)p)0 + ViW > Rn, i = 1,... ,2A: - ^ - l| 

= {p + ?n(p, 0,yi) >0, i = l,...,2k-i-l} , 
we have an obvious upper bound 

I rP+(n(pAyi) 

^ Jo “1“ 

from which, together with (15.5p . we see that 


2k-e-l 

n 

2=1 


dr > <1 


2^^-i/(||(i?n + a(2?n)p)0 + yi||ei) 


2=1 


f(Rnei} 


2k-£-l . . 

< cil^—-—{p + in{p,(),yi)) 


Rn 

, 7(2fc-£-l) 


2=1 ^ 

< C* (max{p, 1}) 

From the argument thus far, for every p > 0, 6 £ and ||yj|| < k(t + s), i = 1,... ,2k — i — 1, 

the expression under the integral sign in ()7.30l) eventually converges to 

i{y e Efpp.e)] /.i(>(0,y), 

while it possesses an upper bound of the form 

OO 

C*(max{p, i})'^-i+T'(2fe-^) ^ l{0 < p < y) 

m=0 

for sufficiently large n. Because of the restriction in e, it is elementary to check that 

^OO 

/ (max{p, V 1{0 < p < e-^dp < oo . 

do m =0 

As a result of the dominated convergence theorem, we have obtained, as n ^ oo, 

r roo p 

Jk ~ a{Rn)RA fiRneiA^ / J{e)d9 dp dy 

JSd-i Jo J(R‘*)2fe-^-i 

„ i{y 6 Efpp.e)] ft!'](0.y) 

= a(R„)Ri-^f{R„e,f=-‘e\{(k - s), 

where the last step follows from the rotation invariance of h.. Hence, we have 

k 

Cov{Gn,K,L{t), Gn,K,L{s)) ~ n^’^-^a{Rn)RV f (Rneif^-^ Me,K,L{t, s), n^oo. 


e=i 


If n/(i?„ei) ^ 0, then the fcth term in the sum is asymptotically dominant, and therefore, statement 
(i) of the theorem is complete. However, the first term becomes dominant when nf{Rnei) —>■ oo. 
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in which case, statement {Hi) is established. In addition, if nf{Rnei) —^ G (0, oo), all the terms 
in the sum grow at the same rate, and this completes statement {ii). □ 

Subsequently, we show the results on finite-dimensional weak convergence of Xji and (X+, X“) 
defined in (I7.28h and (I7.29h , which somewhat parallel those of Proposition 17.31 The reader may 
return to Section [5.21 to recall the definition and properties of the limit (W^, W^). We omit their 
proofs, since the argument in Proposition 17.31 does apply again with minor modifications. 

Proposition 7.6. Assume the conditions of Theorem 15.31 Then, weak convergences (i) — (Hi) in 
the theorem hold in a finite-dimensional sense. Furthermore, the following results also hold in a 
finite-dimensional sense. 

(i) If nf{Rnei) —>■ 0 as n ^ oo, then 

(X+,X-)^(W+,W,~). 

(a) If nf{Rnei) ^ ^ G (0, oo) as n ^ oo, then 

(X+,X-) ^ 

\e=i £=i 

(Hi) If nf{Rnei) —>■ oo as n ^ oo, then 

(X.+x-)^(w+,wr). 

For the same reason as discussed in the preceding subsection, the next proposition can complete 
the proof of Theorem 15.31 

Proposition 7.7. The sequences (X+) and (X~) are both tight in P[0, oo), regardless of the limit 
of nf{Rnei). 

Proof. We only prove the tightness of (X+) but suppress the superscript “-I-” from the functions 
and objects involved during the proof. Proceeding completely in the same manner as Proposition 
we have only to show that there exists B > 0 such that 

(7.32) t-2e{ E E E E hn,t,s{yi) hn,t,s{y2) hn,s,r{y3) hn,s,r{yi) 

yicVn y2CVn ysoVn yiCVn 

X 11 each Ti has at least one common elements 

with at least one of the other three} | < B{t — r)^ 

for allO<r<s<t<L and n > 1. There are three possibilities to be discussed. 

[I] ^ := |Ti n T 2 I G {1,..., k}, f := IT 3 n T 4 I G {1,... , k}, and (Ti U T 2 ) n (T 3 U T 4 ) = 0- 
m£:= |T 2 nT 3 | G {!,..., A:}, f := |TinT 4 | G A:}, and (T 2 U T 3 ) n (Ti U T 4 ) = 0- 

[IE]. Each Ti has at least one common element with at least one of the other three, but neither [I] 
or [n] is true. 

Let B be the left hand side of ()7.32|) . As for case [I], by mimicking the argument in Proposition 
EH we obtain 

B < C*T-^n^^-^-^'a{Rn)^Rl^‘^-^^f{Rnei)^^-^-^'{t - r)^ < C*{t - rf, 

which proves (17.321) . Since we can deal with [E] in an analogous way, we can turn to case [IE]. 
Letting I := Ak — |Ti U T 2 U T 3 U T 4 I G {3,... , 3A:}, the same argument as Proposition 17.41 yields 

B < C*T-V-U{Rn)Ri-^f{Rnei)^'^-\t - rf <C*{t- rf 



which verifies (17.321) . 


□ 
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8. Appendix 

We collect supplemental but important results for the completion of the main theorems. This 
result is known as the Palm theory of Poisson point processes, which is applied a number of times 
throughout the proof. 

Lemma 8.1. (Palm theory for Poisson point processes, [2], Corollary B.2 in [9], see also Theorem 
1.6 in [ 2 T]j Let (Xi) be i.i.d. W^-valued random variables with common density f. Let Vn be 
a Poisson point process on with intensity nf. Let h{y), hi{y), i = 1,2,3,4 he measurable 
bounded functions defined for y S Then, 

e{ e '■(u} = ^ my )}. 

ycVu 

and for every G { 0 ,..., k}, 

2k—£ 

E E hi{yi)h2{y2)i{\yiny2\ = i}} = —^^-^E[hi{yi)h2{y2)i{\yiny2\ = i}}. 

yiCVr,y2CVr, 

Moreover, for every G {0,..., A;} and £ G {0,..., 2A:}, there exists a constant C > 0, which 

depends only on ii, I, and k such that 

E E hi{yi)h2iy2)h3{y3) 

yid.'Pn 3^2C7^7). ysCiT^Ti 

X i{|TinT 2 | =4, =^ 2 , | 3 ^ 2 nT 3 | =4, |3^i u T 2 u T 3 I = 3 fc - ^}| 

= Cn3'=-^]E{/li(Tl)/l2(3^2)h3(T3) 

X i{|TinT 2 | =^ 1 , ITinTsI =^ 2 , |T2nT3| =4, |Ti u T 2 u T 3 I = 3A: - . 

Similarly, for iij G {0,... ,k}, rUp^g^r € {0,.. .,k}, i,j,p,q,r G {1,2,3,4} with i ^ j,p ^ q,p ^ 
r,q r, and i £ {0,... ,3k}, there exists a constant C > 0, which depends only on iij, nip^q^r, 
and k such that 

4e E E E hiiyi)h2iy2)h3{y3)h4iy^) 

y2CPn y^aVn yAdVn 

X l{ |Ti hi I = ^i,j, i, j £ { 1 , 2 , 3,4}, i ^ j , 

|Tp nTg n w| = rup^g^r, P,q,r £ { 1 , 2 , 3 ,4}, p^q,p^q,q^r, 

IT 1 UT 2 UT 3 UT 4 I =4A:-^}| 

= Cn^'=-%{hi(Ti) h2{y2) hsiys) /i4(3^4) 

X l{ |Ti hi 3{; I = ^i,j, i, j £ { 1 , 2 , 3,4}, i ^ j , 

|3^p nTg nwl = mp^g^r, P,q,r £ {1,2,3,4}, p^q,p^q,q^r, 

IT 1 UT 2 UT 3 UT 4 I = 4A:-^}|. 
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